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High temperature expansions for the free energy, the susceptibility and the second correlation moment of the 
classical iV-vector model [also denoted as the O(N) symmetric classical spin Heisenberg model or as the lattice 
O(N) nonlinear sigma model] have been extended to order /3 21 on the simple cubic and the body centered cubic 
lattices, for arbitrary N. The series for the second field derivative of the susceptibility has been extended to order 
f3 17 . An analysis of the newly computed series yields updated estimates of the model's critical parameters in good 
agreement with present renormalization group estimates. 



1. Introduction 

We have extended to order 21 the computa- 
tion of the High Temperature (HT) series for the 
free energy, the susceptibility x(/3,7V), the sec- 
ond correlation moment /j,2(/3,N) and the second 
held derivative of the susceptibility X4,(fl> N) for 
the classical iV-vector model with any TV and on 
all bipartite lattices, e.g. the (hyper) simple cubic 
(sc) and (hyper) body centered cubic (bcc). The 
HT coefficients are written as explicit (rational) 
functions of N and it is also possible to exhibit 
their (polynomial) dependence on the space di- 
mensionality d. 

The model Hamiltonian is 

H{v} = -1/2 v{x)-v{x') 

where v(x) is a A^-component classical unit 
spin. As well known, for N = we get the self- 
avoiding walk model, for N = 1 the Ising 5 = 1/2 
model, for N = 2 the XY model, for N = i the 
classical Heisenberg model, for A~ — > oo the (ex- 
actly solvable) spherical model. We have used 
the vertex-renormalized Linked Cluster Expan- 
sion (LCE) method. Due to lack of space we ad- 
dress the reader to Refs. f|-|} for an extensive bib- 
liography on the LCE technique and on the HT 
series published and analyzed before our work. 
Here we shall only quote the paper by M. Liischer 
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and P.Weisz Q who computed HT series by LCE 
through /3 14 for x, fi2 and \4 on the sc lattices 
in d = 2,3,4. Our work has been made possible 
by an extensive redesign of the algorithms pro- 
posed in Ref. j5| in order to reduce drastically 
the growth of the computational weight with the 
order of the expansion. It is well known that 
the complexity of the HT expansions grows ex- 
ponentially with the order so that obtaining one 
more series coefficient requires an effort substan- 
tially larger than for all previous ones. Let us 
now sketch the calculation of x through order f3 L . 
First we have to list all topologically inequivalent 
multigraphs of the class S2(L) defined as follows: 
i) they have 2 external lines and no more than 
L internal lines; ii) they are connected, 1-vertex, 
1-line irreducible, and contain no odd loops; iii) 
they have only even degree vertices. The Feyn- 
man rules are: for each graph G £ S2 (L) we have 
to compute 

F(G ] N ]d] p) = I J^^4^^ 
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where S(G) is the symmetry number of the graph 
G (an integer number depending only on the 
structure of G); C(G,N) is the group coefficient 
of G, related to the O(N) symmetry group of 
the model (a polynomial in N with integer co- 
efficients); I{G,d) is an integer number counting 
the unrestricted lattice per-site embeddings of G 
(it depends on G and the lattice structure and 
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dimensionality d); W(G, N,d, f3) is associated to 
the vertices of G: it depends on the structure 
of the spin interaction and accounts (" renormal- 
ization") for all vertex insertions in the graphs 
G S2(L) (it is a polynomial in /3 and d and a 
rational function of N). Then we sum over all 
F{G-N-d-(3) 
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and finally 
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It is now clear that in order to get long se- 
ries we need fast and efficient combinatorial al- 
gorithms for producing and comparing graphs, 
since the sum extends to all topologically inequiv- 
alent graphs in ^(i). But generation algorithms 
inevitably produce duplicate graphs which must 
be discarded and identification algorithms have 
a computational weight exponentially increasing 
with the number of vertices! Fast algorithms are 
also needed to compute the symmetry number 
S(G) of each graph: this is also a task of expo- 
nential complexity. Moreover the computational 
load of the algorithm for the embedding number 
I(G, d) grows fastly with d. 

The intricacies and the size of this calcula- 
tion can be gauged from the number of graphs 

2 • 10 7 , which enter into the computation of x 
through order L = 21. This figure should be com- 
pared with the corresponding one w 7 • 10 3 of the 
0(/3 14 ) computation in Ref.||. However we are 
still far from our present computational limits! 

It is also worth to stress that the LCE gives im- 
mediate access also to series for the more general 
model described by the partition function: 

(i,5) 

where tp^ is a iV-component vector. By a proper 
choice of the single spin measure dfj,(tpf) we ob- 
tain the HT series for a variety of models includ- 
ing the general scalar isovector P(</? 2 ) lattice field 
theory, the general spin S Ising model, the Blume- 
Capel model, the double Gaussian model..., and 



this enables us to study various representatives 
of each universality class. The expected returns 
of our laborious enterprise include more accurate 
universality tests, comparisons with estimates of 
critical exponents and of critical amplitude ratios 
obtained within the Renormalization Group (RG) 
approach H either by the Fisher- Wilson pertur- 
bative e-expansion (e = 4 — d) of the continuum 
(if 2 ) 2 model, or by the Parisi coupling constant 
expansion of the same model in fixed dimension 
d = 3, or else by other numerical methods. We 
also expect that the new feature of our calcula- 
tion, namely the explicit dependence on N and d 
of our HT coefficients, can provide further insight 
into the properties of 1/N and e expansions. 

This project has been carried on by an IBM 
Rise 6000/580 w-station with 128 Mbytes mem- 
ory and 1.5 Gbytes disk storage. Typical cpu 
times are a few hours. 

2. Series Analysis 

Our HT series on the square lattice have been 
tabulated and examined in Refs.^Q, while a 
study of the 3d case has been partly presented 
in Refs.0,§|. In the 2d case, our results can 
be summarized as follows: the HT series ex- 
hibit in the range N < 2 a qualitative behavior 
which is sharply different from that for N > 2. 
For N — 2 the critical properties of the model 
are completely consistent with the current ideas 
on the Kosterlitz-Thouless phase transition. For 
N > 2 they appear to be consistent with the con- 
ventional asymptotic freedom expectations from 
the perturbative RG analysis and with exact (al- 
though not completely rigorous) results from the 
Bethe Ansatz. Similar conclusions from an inde- 
pendent computation of 2d HT expansions have 
been reached in Ref.0. In the 3d case, we have 
a more conventional scenario: for all TV a second 
order phase transition occurs at nonzero temper- 
ature. In particular, if we set r = 1 — /3//3 c , in 
the vicinity of the critical point x anc > £ behave 
as 

x(N,0)~A x (N)T-~<W(l + a x (N)T W + 



a' x (N) 



r 26(N) 



+ ■■ + e x (N)r + ..) 
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and 



aN,P)^A i (N)r-^ N \l + a^Ny, 



.e(N)_ 



4(7V)r 2£, W + .. + e ? (7V)r + ..) 

when r | 0. The critical exponents j(N), v(N) 
and the scaling correction exponent 6(N) are uni- 
versal (for each N). On the contrary the crit- 
ical amplitudes A X (N), a x (N), .., A e (N), o^(JV).. 
are nonuniversal, but amplitude ratios like 
a^(N)/a x (N) etc. are universal and interesting 
since they are not yet accurately known. 

The numerical problem of determining simul- 
taneously the critical parameters f3 c (N),j(N), 
v{N) etc. as well as the leading correction expo- 
nent 0(N) is a difficult job: it amounts to an in- 
trinsically unstable double exponential fit. How- 
ever it must be faced, because without proper al- 
lowance at least for the main corrections to scal- 
ing no improvement of the accuracy in the deter- 
mination of critical parameters can be warranted 
even from extended series. We have therefore per- 
formed two kinds of analysis of the series by the 
differential approximant || method: 

i) an unbiased analysis, where we do not as- 
sume to know some critical parameters when try- 
ing to determine the remaining ones; 

ii) a biased analysis, where we assume that the 
exponents of the confluent corrections to scaling, 
take the values indicated by the fixed dimension 
RG perturbative calculations in Ref. |J for N < 3 
and in Ref.@ for N > 3. 

From the analysis of \ and £, we obtain in 
this way numerical estimates of critical param- 
eters which are shown in two extensive tables of 
Ref. jlj, but cannot be reported here for lack of 
space. We can however summarize the qualita- 
tive conclusions. As a first general result, both 
analyses confirm the traditional expectation that 
the bcc series have much better convergence prop- 
erties than the sc series. The results of the unbi- 
ased series analysis agree well with the exponent 
estimates from the fifth order e-expansion over 
the range < N < 3 on which they are available. 
Analogously, we observe an even closer agreement 
with the results from the fixed dimension coupling 
constant six loop expansion |6|,[l0| over the whole 
range of values of N, and, for < N < 3, with 



the seven loop results . However for N > 4, the 
presence of residual trends in the extrapolations 
of the series estimates for the exponents which 
use an increasing number of HT coefficients, seri- 
ously questions the accuracy of the unbiased anal- 
ysis suggesting that the confluent corrections are 
large and inadequately accounted for. Therefore 
a biased analysis is necessary. 

From the biased analysis we also obtain in- 
teresting and, in some cases, accurate estimates 
of various critical amplitudes which confirm that 
scaling corrections should not be neglected. For 
N < 3, our biased exponent estimates essentially 
agree with the unbiased ones, but they have a 
greater accuracy. For N > 4, due to the large 
confluent corrections, the biased estimates of the 
critical exponents differ up to a few percents both 
from the unbiased ones and from the fixed dimen- 
sion six loop perturbation results, suggesting that 
both a seven loop computation and a more accu- 
rate evaluation of the renormalized coupling may 
be needed. 
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